Two-Particle Microrheology of quasi-2D Viscous Systems 
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We study the spatially correlated motions of colloidal particles in a quasi-2D system (Human 
Serum Albumin (HSA) protein molecules at an air-water interface) for different surface viscosities 
r]s- We observe a transition in the behavior of the correlated motion, from 2-D interface dominated 
at high rjs to bulk fluid-dependent at low rjs- The correlated motions can be scaled onto a master 
curve which captures the features of this transition. This master curve also characterizes the spatial 
dependence of the flow field of a viscous interface in response to a force. The scale factors used 
for the master curve allow for the calculation of the surface viscosity rjs that can be compared to 
one-particle measurements. 

PACS numbers: 83.10.Mj, 87.68.-fz, 87.16.Dg 
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Diffusion in three dimensions has been well understood 
since 1905, when two authors showed that the motion of 
particles suspended in a fluid is related to the fluid's vis- 
cosity [l|, Q • This observation has been generalized in a 
technique called microrheology, which measures the ther- 
mal motion of tracer particles introduced in a viscoelastic 
material. From the motions of the particles, the mate- 
rial dependent properties can be determined, such as the 
elastic modulus, G'{lo), and the viscous modulus, G"{llj) 
Q. This has been apphed to measure the viscoelasticity 
of bulk materials such as polymer solutions Q] , biomate- 
rials I3] and hydrogels p- A closely related question is 
the motion of tracer particles in a two-dimensional sys- 
tem such as lipid molecules at an air- water interface JJ or 
lipid rafts in cell membranes [g- For example, in a purely 
viscous 2-D system, one might imagine that the diffusive 
properties are related to the two-dimensional viscosity, 
and that by following the motion of tracer particles one 
could determine this viscosity. However, in most cases 
of practical interest, a strictly two-dimensional surface is 
an idealization and in reality the surface is adjacent to 
three-dimensional fluid reservoirs. For example, recent 
experiments study diffusion in biological systems such as 
cell membranes |8|, |9| which are surrounding a 3-D cell 
and immersed in a 3-D fluid. This coupling modifles the 
behavior of tracer par ticles and makes interpretation of 
the results trickier 



par ticles an 



Furthermore, in many cases in 3D, tracers are known 
to modify the structure of the medium in their vicinity, 
leading to erroneous measurements of rheological quan- 
tities [131 • Another possibility is that pre-existing in- 
homogeneities such as pores in an otherwise rigid mate- 
rial can entrain the tracers, resulting in measurements 
that underestimate the bulk viscoelasticity of the mate- 
rial in question jj]. To overcome these difSculties, a new 
method known as two-particle microrheology has been 
established J^J , which looks at the cross-correlated ther- 
mal motions of pairs of particles. The correlated motion 
of two beads is driven by long-wavelength modes in the 
system, and is therefore independent of the local envi- 
ronment of the tracers. While two-particle microrheol- 
ogy has been applied to 3-D systems [ij, |lj| where the 



strain fleld decays as 1/R, a formalism for 2-D interfa- 
cial systems coupled to viscous bulk fluids has not been 
experimentally determined to date. This is largely due 
to the non-trivial nature of the flow fleld created by the 
thermal motion of particles at an interface. An under- 
standing of this flow fleld is critical towards determining 
the true microrheological behavior of an interface 0)13131 ■ 

In this Letter, we look at the correlated motions of par- 
ticles embedded at an air-water interface in the presence 
of human serum albumin (HSA) protein molecules. This 
is done as a function of particle separation R and lag 
times T, for correlated motion along the line joining the 
centers of particles, Drr{R,T), and motion perpendicu- 
lar, D0g{R, t). The interface is purely viscous, with Drr, 
Dee ~ ■''• The correlations show a transition as a func- 
tion of the surface viscosity rjs ; at high r]s their behavior 
is 2-D dominated; at intermediate rjs they show crossover 
behavior; and at low rjs their behavior is strongly influ- 
enced by the 3-D fluid reservoirs (~ 1/R and 1/R^ re- 
spectively). The correlated motion of the tracer particles 
for different surface viscosities can be scaled onto a single 
master curve, which agrees with theoretical predictions 
ll5J. The surface viscosity rjs determined from the scal- 
ing parameters of the master curve agrees well with r]s 
from one-particle measurements, demonstrating the ho- 
mogeneity of HSA at an interface. 

We use aqueous solutions of HSA over a narrow range 
of bulk concentrations (c = 0.03 — 0.045 mg/ml) to ob- 
tain our interface. At these bulk concentrations, HSA 
molecules diffuse to the air-water interface to form a 
thin monolayer of size ~ 3 nm [l6j, thereby creating 
a surface shear viscosity. The surface concentration of 
HSA slowly increases over time [13, and so the sur- 
face shear viscosity rjs can be varied over a wide range. 
The viscosity of the bulk solution is negligibly differ- 
ent from the viscosity of water, and is assumed to be 
77 = 1.0 mPa-s for all our experiments, while the vis- 
cosity of air is considered to be negligible. Micron- 
sized polystyrene beads (Interfacial Dynamics Corpora- 
tion, carboxyl-modifled, radius a = 0.9/im) dispersed in 
an aqueous solution with 20 % isopropyl alcohol, are 
spread at the interface with a syringe needle to act as 



tracer particles. The particles are imaged by bright-field 
microscopy with a 20x objective, NA=0.45, at a spatial 
resolution of 606 nm/pixel and a frame rate of 30 Hz. 
For each sample, short movies of 200 frames are recorded 
with a CCD camera(Integrated Design Tools, X-Stream 
Vision XS-3) that has 1240 x 1024 pixel resolution, with 
hundreds of particles lying within the field of view. The 
positions of the tracers for every frame are determined by 
particle tracking. Of potential concern is the presence of 
interactions between particles [Tg, whether due to elec- 
trostatic forces, capillary forces, or other origin. For all 
experiments, we measure the pair correlation function 
g{R) and find no structure for the separations R consid- 
ered in our results. Also, our experiments are conducted 
in the dilute particle limit, and a control experiment for 
one of the surface protein concentrations verifies that our 
results are unchanged when the particle concentration is 
varied by a factor of 4. This ensures that long-range par- 
ticle interactions are not present and do not affect our 
conclusions. 

From the particle positions, we determine the vector 
displacements of the tracers Ar^ (t, r) = Va (t + r)— r^ (t) 
where t is the absolute time and r is the lag time. Care 
is taken to eliminate any global drift of the sample from 
these vector displacements. We can then calculate the 
ensemble-averaged cross-correlated particle motions [l^ 

D^p{r,T) ^ {Arlit,T)Ar^^{t,T)6[r - R'Ht)]),^,., (1) 

where i,j are particle indices, a and (3 represent different 
co-ordinates and i?*-' is the distance between particles i 
and j. In particular, we focus on the diagonal elements of 
this tensor product: Drr, which measures the correlated 
motion along the line joining the centers of particles, and 
Dgg , which measures the correlated motion perpendicular 
to this line (the off-diagonal elements are assumed to be 
uncorrelated, and hence 0). In addition, we also calculate 
the one-particle MSD, (Ar^(T)), from the self terms {i = 
j) in the above expression. In 2-D interfacial systems 
coupled to a bulk fluid reservoir, this one-particle MSD 
is related to the surface viscosity r/s by a modified Stokes- 
Einstein relation p^lTlLll^ 



(Ar2(r))=4D;i 



(2) 



where D'^ = D^ [ln{2ris/'na) - 7b + 0{ria/ris)] , with D^ = 
kBT/4:TTr]s and je — 0.577 being Euler's constant. This 
equation, derived by Saffman for rjs/rja ^ 1 [l^, has 
been modified by Hughes [l^ to work in the limit where 
ris/rja ^ 1. Equation 2 has been used to describe a variety 
of homogeneous systems ^20] and we use the measured 
values of (Ar^(T)) to solve it for the one-particle surface 
viscosity, 77s,ip (refer inset of Fig. 1 for example). 

However, as noted above, one-particle measurements 
may be inaccurate if the system is heterogeneous, so we 
perform two-particle measurements to check this, and 
also to further probe the spatial flow field around the 
particles. In Fig. 1(a), we show D^r as a function of R 
for different lag times r, for a sample with r]s,ip = 340 




FIG. 1: (a) and (b) Two point correlation functions Drr(R, t) 
and Dee{R, r) for a sample with »?s,ip = 340 nPa-s-m, with lag 
times r = 0.1, 0.2, 0.4 and 0.8 s. The evenly spaced correlation 
functions imply that Drr, Dee ~ t. Inset: One-particle MSD 
for the sample (diamonds), where the straight line is a fit 
giving 77s, ip = 340 nPa-s-m, using Eqn. 2. 



nPa-s-m. The motion of a tracer particle creates a flow 
field that affects the motion of other particles. This flow 
field decays as we move further out from the particle; 
hence the correlated motions decay as a function of par- 
ticle separation for all lag times. From the one-particle 
measurements, we observe that HSA monolayers at an 
interface are entirely viscous at these surface concentra- 
tions, and therefore we also observe that Drr ^ t. This 
is illustrated in Fig. 1(a) where all the Drr are spaced 
evenly on a log-plot for lag times r = 0.1,0.2,0.4 and 
0.8s. A similar linear relationship exists for Dgg, shown 
in Fig. 1(b). The linear scaling of the correlation func- 
tions enables the estimation of r-independent quantities 
{Drr/T)r and {Dgg/T)r, which depend only on R, and 
have units of a diffusion coefflcient. 

Changing the surface shear viscosity 77s , has a substan- 
tial effect on both Drr and Dgg. At high rjs, Fig. 2(a), 
both (Drr/T) and (Dgg/r) are nearly equal and constant 
over the length scales 10 < i? < 150/zm. At these sur- 
face viscosities, the behavior can be considered to be 2-D 
dominated, where the bulk fluid reservoirs have minimal 
influence. As rjs is decreased, Fig. 2(b), we see curvature 
in these functions and deviation between their values, 
over the same range of length scales. At very low rjs, 
Fig. 2(c), this deviation is even more pronounced, with 
(Drr/r) ~ l/R and {Dgg/r) ^ l/R^. At these viscosi- 
ties, bulk fluid effects begin to dominate, although the 
behavior is still 2-D as the protein molecules are conflned 
to an interface. For all surface viscosities, the behavior of 



the correlation functions is markedly different from what 
is seen in bulk 3-D systems (where D^r, Dee ~ 1/-R)> s^d 
is sensitively dependent on ry^. Despite the differences 
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FIG. 2: Correlation functions (Drr/T") (solid symbols) and 
{Dee /t) {open symbols) as a function of particle separation 
R, for decreasing surface viscosities 77s . The samples are a) 2- 
D dominated: 77s, ip — 340 nPa-s-m (diamonds), b) Crossover: 
rjs,ip = 72 nPa-s-m (circles), c) 3-D dependent: ris,ip = 2f.3 
nPa-s-m (squares). 



where the H^, are Struve functions, and the Y^, are Bessel 
functions of the second kind. More importantly, up to 
a scale factor, Eqn. 3, and therefore the master curve, 
also characterizes the spatial dependence of the flow field 
at an interface in response to a perturbation, such as 
thermal motion of tracers. To our knowledge, this is the 
first experimental mapping of this flow field over such 
a wide range of length scales. The inset to Fig. 3 de- 
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FIG. 3: Master curve of scaled variables _Drr (solid symbols) 
and Dee (open symbols) as a function of reduced particle sep- 
aration p. Symbols are same as in Fig. 2, with two addi- 
tional data sets (»?s,ip = 2.1 nPa-s-m; stars and ?7s,ip = 1275 
nPa-s-m; triangles). Solid lines represent theoretical fits to 
the data. Inset: One- and two-particle surface viscosities for 
the five samples shown in the master curve. The straight 
line has a slope of 1, indicating an equality between the two 
viscosities. 



in the behavior of (Drr/T) and (Dee/r) at low and high 
surface viscosities, all the data can be scaled onto a single 
master curve. We define dimensionless correlation func- 
tions Drr,ee = {Drrfie/T) /2Ds and a reduced separation 
(3 = R/L, where L — rjs/rj. Both scale factors L and 
Dg depend only on r]s; however, we allow them to vary 
independently to obtain two independent measures of 77s 
(77s ^iid 77" respectively, in other words, /? = Rrj/r]'^, and 
Ds = kBT/Anrj'^). Fig. 3 shows the scaled variables Dm 
Dee plotted against the scaled separation /3. All the data 
sets fall on a single master curve, that spans nearly 4 or- 
ders of magnitude. This master curve has the character- 
istics of the individual data sets: at small /? (/3 ^ 1), the 
curves are nearly logarithmic, at intermediate /?(/?« 1), 
they show crossover behavior, while at large (3 {(3 ^ 1), 
they show behavior that asymptotically approaches 1/ R 
and 1/R^ for D^r and Dee respectively. The solid lines 
are fits obtained from theoretical calculations of the re- 
sponse of an interface to an in-plane point force [13, [l3 ■> 
and have the form 

Drr = [f i/l(/3) -j,- f (ro(/3) + Y^im 

f(ro(/3)-F2(/?))l3) 



Dee = [nHoiP) - fMiiP) 



2 



scribes the scale factors used to create the master curve. 
For all the samples, the two independent estimates of 
the two-particle viscosity, 77^ and 77", are within 15% of 
each other; we therefore plot their average, 77^,2^, against 
the one particle viscosity 77^,1^. For most of the samples, 
the one- and two-particle measurements agree reasonably 
well with each other. However, at the highest viscosity, 
and therefore, the highest surface concentration of HSA, 
the two measurements deviate well beyond our experi- 
mental error. We believe this is a consequence of hetero- 
geneity in the system; one possibility could be the for- 
mation of condensed phases at the interface at such high 
surface concentrations of HSA. While this would explain 
the underestimation of the true surface viscosity by the 
one-particle measurement, such an assumption deserves 
further study. 

Conventional microrheology uses the one-particle 
MSD, (Ar2(r)), to determine G'{lu) and G"(w) [I; these 
quantities are directly compared to bulk rheological mea- 
surements. Since the one-particle measurements are in- 
herently local in nature, a more accurate approach is to 
determine the two-particle MSD |l3|, (Ar^(T))D- Anal- 
ogous to [13 , we calculate this quantity for our 2-D vis- 
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FIG. 4: One-particle (solid symbols) and two-particle (open 
symbols) MSDs for four of the five samples shown in Fig. 3. 
Inset: radial scaling of Drr for the different samples at a lag 
time r = 0.53s. 



cous systems by extrapolating the long wavelength fluc- 
tuations of the medium to the bead size, which gives 



(Ar2(r))D 






)R[ln{2rjs/va) - Je + 0(r/a/ry,)](4) 



This expression can readily be generalized for viscoelastic 
interfaces, i.e, with a non-zero G'{uj); however, such mea- 
surements are beyond the scope of this study. In practice, 
we confirm that Drr{R,T)/D{p) is nearly constant over 
the length scales studied, 9 < i? < 100/im, shown in the 
inset to Fig. 4 for a specific lag time (r — 0.53s). The 
averaged quantity {Drr{R,T) / D{fi)) r is then calculated 
for all lag times to obtain {/S.r^ [t)) d . From Fig. 4, we see 
that for all the samples, {/S.r'^{T))o is purely diffusive, as 
expected for a viscous system. At short lag times, the 



turnover in the one-particle MSDs, caused by resolution- 
limited noise, is significantly reduced in the two-particle 
MSDs. Finally, at long lag times, {/\r'^{T))o is equal to 
(Ar^(T)) within experimental error, except for the high- 
est viscosity (as expected, since "qa.ip deviates from 7]s,2p 
for this sample). These observations provide conclusive 
evidence of the accuracy of two-particle measurements 
over local probes of the rheology. 

The verification of two-point microrheological tech- 
niques for a quasi-2D systems has applications for the 
study of inhomogeneous materials at an interface. Any 
significant variation between {/S.r^ (t)) £, and {/S.r^ir)) in- 
dicates the presence of heterogeneities, and the estima- 
tion of rheological quantities from the motion of tracer 
particles can be modified to reflect this. Future work 
will involve studying lipid molecules at an interface un- 
der compression/expansion, which creates domains such 
as liquid expanded or liquid condensed phases. We ex- 
pect that our two-particle measurements will provide ac- 
curate measurements of the surface viscosity even with 
the presence of these heterogeneities. Other possible ar- 
eas of future research are biological interfaces such as 
cell membranes, with the diffusing entities being protein 
aggregates or lipid rafts. However, two-particle microrhc- 
ology can only probe heterogeneities of the order of the 
bead size and above, which may not always be appli- 
cable for lipid systems with molecules in the nanometer 
scale. It should also be pointed out that for cases where 
?7s is very small, the correlation functions [Dgg in partic- 
ular) die out rapidly, making an estimation of ?7s,2p from 
the two-particle measurements difficult. However, one- 
particle microrheology is equally inaccurate in this limit, 
as the bulk viscosity will dominate over surface effects, 
making such measurements challenging. Nonetheless, it 
would be extremely interesting to test the limits of the 
scaling behavior in this regime. 
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